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^ 1 INTRODUCTION 



ABSTRACT 

A new component of the cosmic medium, a light scalar field or "quintessence ", has 
been proposed recently to explain cosmic acceleration with a dynamical cosmological 
constant. Such a field is expected to be coupled explicitely to ordinary matter, unless 
some unknown symmetry prevents it. I investigate the cosmological consequences of 
such a coupled quintessence (CQ) model, assuming an exponential potential and a 
linear coupling. This model is conformally equivalent to Brans-Dicke Lagrangians with 
power-law potential. I evaluate the density perturbations on the cosmic microwave 
background and on the galaxy distribution at the present and derive bounds on the 
coupling constant from the comparison with observational data. 

A novel feature of CQ is that during the matter dominated era the scalar field has a 
finite and almost constant energy density. This epoch, denoted as 0MDE, is responsible 
of several differences with respect to uncoupled quintessence: the multipole spectrum 
of the microwave background is tilted at large angles, the acoustic peaks are shifted, 
their amplitude is changed, and the present 8Mpc//i density variance is diminished. 
The present data constrain the dimensionless coupling constant to < 0.1. 



X 



The recent evidence in favour of an accelerated cosmic expansion (Perlmutter et al. 1999; Riess et al. 1999) has prompted the 
theorists to hypothesize components of the cosmic medium additional to the ordinary matter and radiation, whose equation 
of state is unable to provide the required kinematics. In a flat Universe, the dark energy of such a component should provide 
roughly 70% of the cosmic density, and should possess an effective equation of state 



p = (™ - l)p, 

with the present value (Turner & White 1997, Perlmutter et al. 1999, Wang et al. 1999) 
w G (0,0.6). 



(1) 



(2) 



The most obvious candidate, a cosmological constant, which provides w = 0, has unappealing features: its value would 
be one hundred orders of magnitude smaller than dimensionally expected; upper limits from lensing effects barely allows for 
a f^A = 0.7 (Kochanek 1995), as would be necessary to reconcile the amount of matter in clusters with the ffatness suggested 
by inflation. The next simplest possibility is perhaps to include in the cosmic fluid a light scalar field. In fact, if the field is 
fight enough to vary slowly during a Hubble time, its potential energy can drive an accelerated expansion, just like during 
inflation. The varying field equation of state can then be tuned to lie in the observed range: if this is the case, then the scalar 
field is sometimes denoted in the literature as "quintessence". The scalar field density fraction can be made to decrease 
rapidly in the past, so as to pass easily the lensing constraints, and to avoid discrepancies in the primordial nucleosynthesis 
abundances. 

Beside the acceleration argument, a light scalar field is interesting on its own. First, it is predicted by many fundamental 
theories (string theory, pseudo-Nambu-Goldstone model, Brans-Dicke theory etc.), so that it is natural to look at its cosmo- 
logical consequences (Ratra & Peebles 1988, Wetterich 1995, Frieman et al. 1995, Ferreira & Joyce 1998). Second, the presence 
of a scalar field may fix the standard CDM spectrum (Zlatev et al. 1998, Caldwell et al. 1998, Perrotta & Baccigalupi 1999, 
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Viana & Liddle 1998). Finally, even a small amount of scalar field density may give a detectable contribution to the standard 
CDM scenario, similar to what one has in the MDM model (Ferreira & Joyce 1998). 

A scalar field, however, is expected to couple explicitly (that is, beyond the gravitational coupling) to ordinary matter, 
with a strength comparable to gravity, as put into evidence by Carroll (1998), unless some special symmetry prevents or 
suppresses the coupling. Such a strong coupling would render the scalar field interaction as strong as gravity, and would 
therefore have been already detected. However, a residual coupling still below detection cannot be excluded; moreover, if the 
coupling to baryons is different from the coupling to dark matter, as proposed by Damour et al. (1990), then even a strong 
coupling is indeed possible. Exactly the same arguments hold if one supposes the quintessence field to be coupled to gravity, 
rather than to matter, as investigated by Uzan (1999), Chen & Kamionkowsky (1999) and Baccigalupi et al. (1999). Indeed, 
the two models, although physically different, are related mathematically by a conformal transformation (Wetterich 1995, 
Amendola 1999a). 

The non-minimal coupling of the quintessence field to ordinary matter is therefore worth investigating, especially because 
the wealth of high-precision data that is near to come allows the intriguing possibility of detecting the coupling on the 
microwave background and on the present galaxy distribution. In a previous paper (Amendola 1999b, hereinafter Paper I) 1 
showed that a scalar field with an exponential potential (Wetterich 1988, Ratra & Peebles 1988) and an explicit coupling to 
matter may behave as a kind of hot dark matter component, as was first shown by Ferreira & Joyce (1998) for zero coupling. 
I showed that the CMB spectrum of the model presents acoustic peaks displaced from their location without coupling, and 
that the galaxy power spectrum also bends in agreement to real data. In that case, the field density amounts to at most 20% 
of the critical density, and the expansion is not accelerated. The interesting feature was that the universe has always been in 
an attractor solution, independently of the initial conditions. 

In this paper I focus instead on accelerated solutions. I explore first the general phase space of a homogeneous quintessence 
model with the same exponential potential and coupling to matter as in Paper I; I will refer to this model as coupled 
quintessence, or CQ. Once the phase-space attractors have been identified, two distinct solutions are selected that allow an 
accelerated epoch, and the density fluctuations on these trajectories are studied by the use of a purposedly modified version of 
the code CMBFAST by Seljak & Zaldarriaga (1995). The linear perturbations in the uncoupled case has been already studied 
by Viana & Liddle (1998) and Caldwell et al. (1998). As it will be shown, the coupling introduces several qualitatively new 
features. 



2 COUPLED SCALAR FIELD MODEL 

Consider two components, a scalar field (j) and ordinary matter (e.g., baryons plus CDM) described by the energy-momentum 
tensors T^,ji^^) and r,,^(m,), respectively. General covariance requires the conservation of their sum, so that it is possible to 
consider a coupling such that, for instance, 

^^(m);^ = ~CT(,^)Cj>;^. (3) 

Such a coupling arises for instance in string theory, or after a conformal transformation of Brans-Dicke theory (Wetterich 
1995, Amendola 1999a). It has also been proposed to explain 'fifth- force' experiments, since it corresponds to a new interaction 
that can compete with gravity and be material-dependent. A coupling that violates general covariance is instead studied in 
Barrow & Magueijo (1999). 

The specific coupling (^ is only one of the possible forms. Non- linear couplings as Cr(„)_F((^)(^;^ or more complicate 
functions are also possible. Also, one can think of different couplings to different matter species, for instance coupling the 
scalar field only to dark matter and not to baryons, as proposed by Damour et al. (1990) and Casas et al. (1992). Notice that 
the coupling to radiation (subscript 7) vanishes, since Tj.y) = 0. Here I restrict myself to the simplest possibility (^), which 
is also the same investigated earlier by Wetterich (1995) and is the kind of coupling that arises from Brans-Dicke models. In 
fact, a field with a gravity-coupling term ^^(p'^R in the Lagrangian acquires, after conformal transformation, a coupling to 
matter of the form (^. In the limit of small positive coupling one obtains 

C = ^TC- (4) 

where = SnMp^ and Mp is the Planck mass. Moreover, if the Brans-Dicke Lagrangian contains a power-law potential 
V{<l)) ~ , then the transformed field acquires an exponential potential that, for small positive ^, can be written as 



V{(f)') ~ exp (n — 4) Ky^^^' 



(5) 



(see e.g. Futamase & Maeda 1989, Amendola et al. 1993, Amendola 1999a). The CQ model with a linear coupling and an 
exponential potential that is studied here is therefore conformally equivalent to a large class of Brans-Dicke Lagrangians. 
There are several constraints on the coupling constant C along with constraints on the mass of the scalar field particles. 
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reviewed by Ellis et al. (1989) and Damour (1996). The constraints arise from a variety of observations, ranging from Cavendish- 
type experiments, to primordial nucleosynthesis bounds, to stellar structure, etc. Most of them, however, apply only if the 
scalar field couples universally to all matter, which is not necessarily the case. The most stringent bound for a universal 
coupling, quoted by Wetterich (1995), amounts to 

\C\<0.1Mp\ (6) 

If the coupling to dark matter is different from the baryon coupling, then the constraints on the former rela^KCS considerably, 
and becomes (Damour et al. 1990) 

\C\ < 5M-\ (7) 

It is also to notice that these constraints are local both in space and time, and could be easily escaped by a time-dependent 
coupling constant. In the following, therefore, C is left as a free parameter. 

The constraints from nucleosynthesis refer to the energy density in the scalar component. This has to be small enough 
not to perturb element production, so that at the epoch of nucleosynthesis (Wetterich 1995, Sarkar 1996, Ferreira & Joyce 
1998) 

n^iTns) < 0.1 - 0.2. (8) 
This bound is satisfied by all the interesting models. 



3 BACKGROUND 

Here I derive the background equations in the flat conformal FRW metric 

ds^ = a^{—dT^ + Sijdx'^dx^). 
The CQ scalar field equation is 

4> + 2H4> + a''U,4,^Cp,„a^, (9) 
where H = a/a , and I adopt the exponential potential 

U{^)=Ae\fi'""^. (10) 
The matter (subscript m) and the radiation (subscript 7) equations are 

pm + "iHpm = -Cpm4> (H) 

p^ + AHp^ = 0. (12) 
Denoting with tq the conformal time today, let us put 

O(ro) = 1, Pm(To) = -5 f2m = PmO, Pf{To)=P',0, 0(ro)=0O. (13) 

OTT 

Without loss of generality, the scalar field can be rescaled by a constant quantity, by a suitable redefinition of the potential 
constant A. We put then = 0. This gives 

-3 -Cd> 

Pm = Pmoa e ^, 

Pt = p-,oa^^. (14) 
The Friedman-Einstein equation can be written as 

6 \ a a'' 2 ) 

The dynamics of the CQ model has been analysed in Amendola (1999a) in the regime in which either matter or radiation 
dominates. Here I generalize the analysis to the case in which both matter and radiation are present. As we will see, this 
introduces some interesting new features. Generalizing Copeland et al. (1997) the following variables are introduced: 

along with the independent variable a = log o. Notice that , and give the fraction of total energy density carried by 
the field kinetic energy, the field potential energy, and the radiation, respectively, that is f?,^ = x^ and Clj = z'^. Clearly, 
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the matter energy density fraction is the complement to unity of a;^ + j/^ + 2^. We can rewrite the Eq. (^) and ( |l5|) as 



^ 1 \ 2 I 2 2 



2 



fj.xy + y [2 + ^ 



2' = -|(l-3x^+3y^-2^), (17) 
where the prime denotes here d/da and where I introduced the dimensionless constant 

(in Amendola 1999a P was defined as twice the value above). The parameters /3 and jj, are all we need to completely specify 
our model. The constraints quoted in the previous section on C become now /3 < 0.025 for the universal coupling and /3 < 1 
for the dark matter coupling. 

The system ( [l7| ) is invariant under the change of sign of y, z and of a. Since it is also limited by the condition p > to 
the circle +y^ + < 1, we limit the analysis only to the quarter of unitary sphere of positive y, z. The critical points, i.e. 
the points that verify x' = y' — z' — 0, are scaling solutions, on which the scalar field equation of state is 

2x^ 

^0 = — o = const, (19) 
x^ + y' 

the scalar field total energy density is = x^ + y'^, and the scale factor is 

a ~ rT^ = p= (20) 

(t being the time defined by dt — a{T)dT). The effective equation of state for the total cosmic fluid, ptot = {wctt — l)ptot has 
index 

"Weft — 1 + x'^ — y'^ + 2^/3 = 1 + fl~i(w^i — 1) + r2^(ui0 — 1) 

(where w-y = 4/3 is the radiation equation of state). As already noticed, a value < < 0.6 is required by observations, 
while WcS < 2/3 is enough for acceleration. 

The system ( |l7| ) with an exponential potential has up to fifteen critical points, of which only eight can be in the allowed 
region. They are labelled by a letter that reproduces the classification given in Amendola 1999a, and a subscript that denotes 
whether beside the field there is a component of matter (subscript M), radiation (R), both (RM) or neither of the two (in 
which case the energy density is taken up completely by the scalar field; no subscript in this case). The critical points are 
listed in Tab. I, where g(/3, /i) = 4/3^ + 4/3/i + 18. For every value of the parameters //, /3 there is one and only one stable 
critical point (attractor); one or more saddle points can also exist. More details on the phase space dynamics in Wetterich 
(1995) and Amendola (1999a). 



Point 


X 


y 


2 




p 




W0 


a 


M 
3 







1 


3//i' 


2^iV9 


2^iV9 


bR 


M 


ImI 




c 


1/2 


4/3 


4/3 


bM 


3 


(9-9)^/' 





a 




M 


18 




2|m+/3| 








CR 








1 





1/2 


4/3 




CRM 


1 

2/3 







1 


1/2 


4/3 


2 


CM 


¥ 










6 

4/3^+9 


'-^ 9 


2 


d 


-1 








1 


1/3 


2 


2 


e 


+1 








1 


1/3 


2 


2 


Tab. I. Critical points. 
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The regions of existence and stability of the critical points are summarized in Tab. II . In the table the attention is 
restricted to the half plane /i > 0, since there is complete symmetry under simultaneous sign exchange of fi and /3. I defined 

Fig. 1 displays the regions of the parameter space > in which the various points are stable. 



Point Existence 


Stability 


Acceleration 


a fj, < 3 







bji fi > \/6 < /i < —4/? never 



bM \ti + l3\>3/2,ii>iio A* > M+, A* > -4/3 /i < 2/3 



cr 




unstable V/^, (3 


never 


CRM 


m > V3/2 


fi > -4/3 


never 


Cm 


m < 3/2 


\I3\ < VS/2,fi< fio 


never 


d 




unstable V/i, /3 


never 


e 


Vm,/3 


unstable V/i, /? 


never 



Tab. II. Properties of the critical points. 



There are only two critical points that admit accelerated solutions, i.e. solutions that satisfy (|^: the points a and bM- 
They differ in several important aspects, so we study them separately. 
Solutions of type a. 

The attractor a, once reached, brings to zero the matter density. To allow for the observed matter content of the universe, 
we have to select the initial conditions, if they exist, in such a way that the attractor is not yet reached at the present 
time, but the expansion is already accelerated. On the positive side, this attractor is accelerated also for small values of the 
coupling constant. Before discussing in detail the solutions, let us notice that the limit f3 = fj, — corresponds to the ordinary 
cosmological constant. Suppose then we put initially the field at zero kinetic energy (a; = 0). A trajectory acceptable from a 
cosmological point of view should start deep into the radiation era (z ~ 1), then enter a matter dominate era (z « 0), and 
finally fall into the attractor a, the only one still existent for P — /j. = 0, which corresponds to the A-dominated stage. In 
other words, the path of a ordinary A universe would be cr cm a. A similar sequence of critical points characterizes all 
the trajectories discussed in the following. 

In Fig. 2 we show the 3D phase space {x,y,z) of model a, with jS = 0, = 0.1 and (3 — 0.5, /i — 0.1 . As before, a 
cosmological trajectory must start in the radiation era {z ~ 1) and has to provide the correct final conditions for x, y and 
H. Since the scalar field has to start dominating only recently, it is clear that its initial energy density deep in the radiation 
era has to be very small: in the 3D phase space this means that the cosmic solutions start near {x, y, z) = (0, 0, 1) , that is, 
near the unstable critical point cr. The trajectories in Fig. 2 that fall almost vertically from top are examples of such cosmic 
solutions. The attractor of the CQ model is the same as for the uncoupled case, at (a;, y, z) = (—0.033, 0.99, 0), but while in 
the ordinary quintessence case there is a saddle point cm at {x,y,z) = (0,0,0), in CQ this moves to {x,y,z) = (2/3/3,0,0), 
on which = 4/3^9. This appears more clearly from Fig. 3, in which the trend of 0^,0^ and 0,-^ is reported. The path of 
this solution is then cr cm ~^ a,, just as in the pure A model. For /3 > 3/2, actually, the point cm ceases to exist, but such 
high values of (3 are anyway unacceptable. 

It is to be emphasized that the stage of constant and finite = 4/3^/9 is typical of the CQ model, as it is absent in the 
limit of zero coupling. Let us call this the field-matter-dominated era, or (/)MDE . As we will see, this stage is responsible of 
most of the differences with respect to ordinary quintessence. 

The equivalence time in CQ occurs earlier than in uncoupled quintessence: 



P70 

pmO 



(21) 



For the acceptable values of /3, however, this shift has only a minor effect. 

Finally, it is to be noticed that the attractor a is independent of the coupling (3. Then, as the universe at z ~ 5 converges 
toward the point a, the dynamics becomes independent of the coupling. As a consequence, the cosmological probes at z < 5 , 
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like the type la supernovae or the cluster abundance, are not efficient tools for discriminating between ordinary quintessence 
and CQ. 

Solutions of type &m- 

The attractor bM is a solution for which the matter density and the scalar field share a finite and constant portion of the 
cosmic energy. For instance, if we put 

P = 4.02, = 2.68 (22) 

we get Q,M = 0.3 , Q,^, = 0.7 and ~ 0.14, well within the requested range. These values, once reached, remain indefinitely 
constant. The coincidence of similar values of the energy density in the matter and field component is therefore solved, 
independently of the initial conditions. On the negative side, however, these solutions require a large value of the coupling 
constant { P > -^/B/S to obtain < 0.6) and are therefore at risk of running into conflict with constraints on the coupling 
derived from local experiments. The strongest objection to these solutions, however, comes from the simple fact that they 
lack a matter dominated era, as we will show in a moment. 

In Fig. 4 we show the phase space of model bM, assuming /3 = 4.02, /j, = 2.68. As can be seen, the phase space is 
now completely different. The attractor is at {x,y,z) = (—0.22,0.81,0). The trajectory that falls from top is again similar 
to the one used in the perturbation calculations. As can be seen better in Fig. 5, there is a transient near the saddle point 
CRM , here at {x, y, z) ~ (0.12, 0, 0.95). In Fig. 5 the evolution of n„i, and Q.^ in two cases is displayed, one for which the 
present value of Woff equals 0.4, for which [3 and fj, must be chosen as above, and the other for w^s — 0.5, which requires 
P = fjL = 2.37. In both cases we can see the transient crm , characterized by Q,^ — 1/(4/3^) and fijv/ = 1 — 1/(2/3^), followed 
by the decay of the radiation component and the stabilization of the field and matter to their final values. The path of this 
solution is then cr — > crm — > &m, in contrast with the solutions of the type a. As already remarked, such a trajectory is 
possible only for /3 7^ 0. The most conspicuous features here is that the radiation dominates until recently [z ~ 50) and 
the matter dominated era is absent. As is intuitive, such behavior is catastrophic for the growth of the fiuctuations: when a 
fiuctuation mode reenters the horizon, in fact, is suppressed first by the long RDE, and then by the accelerated expansion. 
As a consequence, the present as is unacceptably small, of the order of lO^"* for COBE normalized spectra. Unless some 
mechanism other than gravitational instability powers the fluctuations, the trajectories of type bM are precluded. In the 
following we restrict therefore our attention to solutions of type a. 

This concludes the analysis of the homogeneous solutions of CQ. It is to be stressed that we considered all the possible 
accelerated solutions. As in the next section we will span all the range of /3 and /i that yield cosmologically acceptable solutions 
of type a, we can consider exhausted the analysis of CQ for exponential potentials and linear coupling. 



4 PERTURBATIONS ALONG SOLUTIONS A 

We now proceed to study the evolution of the perturbations in the coupled quintessence theory. The equations of the per- 
turbations in the synchronous gauge have been derived and discussed at length in Paper I and will not be repeated here. In 
that case it was shown that several important features of the perturbation evolution could be derived analytically, since the 
background evolution was always on one of the attractor, and thus particularly simple. The same occurs here, at least in some 
cases. As described in Paper I, all the results presented below have been obtained by modifying the code CMBFAST of Seljak 
& Zaldarriaga (1996). 

The key fluctuation equation in Paper I was the evolution of the sub-horizon perturbations in the MDE regime, the only 
situation in which the evolution differs from the pure CDM case. For CQ, this regime is actually the <j) MDE introduced above. 
Let us note flrst that along any attractor with x = Xa one has, from Eq. ( [l^ ) 

<^=^^loga. (23) 

Denoting with 5 the fiuctuation in the matter component, we derived in Paper I an equation for sub-horizon modes in MDE 
along an attractor that can be rewritten as follows: 

5 + H{l + 2xal3) S - ^H'^QmS ^1 - = 0. (24) 

As we noticed in the previous section, the solutions a passes through the (^MDE transient attractor cm, before the present 
^-dominated epoch. Therefore, from Tab. I, we can substitute Xa = 2/3/3. It follows that Eq. ( p^ has a growing solution 

5 = Aa"', with m=l-^. (25) 

This shows two important facts: flrst, the fiuctuations in 0MDE are suppressed with respect to the standard CDM behavior 
(m = 1), which also holds for the uncoupled quintessence model, for all values of /3; second, the evolution does not depend on 
the sign of (3. B efore the present time, at 2: ~ 5, the trajectory deviates from the 0MDE solution, and begins to dominate. 
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In Fig. 6 the behavior of three density fluctuation wavelengths, calculated numerically with CMBFAST , is shown as a 
function of the scale factor a. I plot S/a to enhance the differences among the various cases. It is possible to distinguish four 
distinct epochs. Let us first follow the intermediate wavelength of Fig. 6. First, the 100 Mpc//i fluctuation grow as , while 
it is a super-horizon mode in RDE; then, around z ~ 2500 (a ~ .0004) , it reenters the horizon and freezes until the ^MDE 
begins. In uncoupled quintessence (/3 — 0) the fluctuation grows now as a , as usual; in CQ the growth is instead suppressed 
as expected from Eq. (^5|). Finally, around « ~ 5, </) starts to dominate, the universe accelerates, and the fluctuation growth 
is definitely suppressed in all cases. The longer wavelength follows the same phases, except that it reenters directly in the 
<I)MDE, and therefore bypasses the freezing stage. The shorter wavelength starts in the plot already inside the horizon, and 
follows the same (jiMDE growth evolution of the other modes. 

The same ^MDE attractor solution can be used to derive the location of the first acoustic peak on the CMB. In fact, 
this depends essentially on the size of the sound horizon at decoupling (subscript d). We have shown in Paper I that the 
following approximated law governs the size of the sound horizon along an attractor solution 

where ro — 2{ad/'iY^'^ H^^ is the standard sound horizon. Therefore, the multipole location of the first acoustic peak is 



peafe ^ ^lo{l + 2XaP)a-'"^^ (27) 



3 



(28) 



where ~ 2tt/ (roHo) — 200 is the standard peak location, and where the second line specializes to the (/)MDE attractor cm- 
For instance, /3 = 0.1 gives a location Ipeak — 1.06^o which, at least to a first approximation, is in agreement with the numerical 
results below. Notice again that the peak displacement is always toward larger ^s, regardless of the sign of l3, contrary to what 
happened in Paper I along the attractor bM. Similar behavior was found numerically by Chen & Kamionkowsky (1999) and 
Baccigalupi et al. (1999) in Brans-Dicke models. Here, substituting = 3^/2, we find that the peak shifts approximately as 

The solution we use as background in this section is, as anticipated, on its way to the attractor a. The initial conditions 
will be chosen so that at the present time (i.e., when = 3000 Mpc//i) we have Q,m ~ 0.3 , Q,^ — 0.7 and ~ 0. This 
implies that at the present time we should have 

xo ~ 0, 2/0 = 0.837, zo = 0.092 (29) 
independently of /3. We begin by investigating the parameter range: 
M = 0.1, /3 6(0,0.15) 

The initial conditions that produce the requested final values have been obtained by trial and errors. Inserting the background 
solution in the modified CMBFAST code, we obtain the CMB spectra reported in Fig. 7. The other values adopted are 

h = 0.7, fib = 0.04, 71 = 1. (30) 

It can be noticed that the peaks move to larger multipoles, as expected. Their amplitude is generally reduced, due both to 
the growth suppression found above, and to the fact that now the COBE normalization at small £ includes the integrated 
Sachs- Wolfe (ISW) effect, no longer negligible in CQ, and as a consequence the fluctuation amplitude at decoupling is reduced 
(see for instance Hu & White 1996). The ISW is also responsible for the tilt at small multipoles. Models with /3 > 0.15 
are aheady ruled out by CMB observations, while models with /3 < 0.01 are essentially indistinguishable from uncoupled 
quintessence. Values smaller than /3 < 0.06 produce acoustic peaks which are slightly above those for the uncoupled model. 
The bound /9 < 0.15 is already stronger than (Q), valid for the coupling to dark matter; the determination of the spectrum 
with a precision of 5%, within reach of the Boomerang or Maxima experiments will constrain (3 to two decimal digits, better 
than current constraints to the baryon coupling (^. The effect for /3 < 0.03, which satisfies the constraint for the universal 
coupling model, will be distinguishable in the near future. 

The reason for the small increase in the acoustic peak that is observed for |/3| < 0.06 is not entirely clear. Since the 
numerical fluctuation growth follows very closely the analytical prediction of Eq. (^) I believe the rise is not a numerical 
artifact. Notice that for small /? the 0MDE transient attractor starts just past the decoupling epoch, and thus the analytical 
expectations based on the </!> MDE solution are not accurate. 

The effect of changing /i, within the limit /i < \/3 necessary to have acceleration, is minimal, since the trajectories must 
anyway satisfy the same final conditions. The spectrum for /3 = is therefore almost identical to the spectrum of a pure A 
model with f^A = 0.7 . Also, as expected from the analytical expressions, the perturbative results are almost insensitive to 
the sign of /3. The present analysis therefore spans all the possible accelerated trajectories in CQ with Sl,^ = 0.7 that are 
cosmologically acceptable. 
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In Fig. 8 I report the power spectra A^(fc) = fc^P(fc)/(27r^) normalized to COBE, compared to the data compiled and 
corrected for redshift and non-linear distortions by Peacock & Dodds (1994). It can be seen that the slope of the spectrum is in 
rough agreement with what is observed (the largest discrepancy is for the four smallest scale data points, where non-linearity 
and redshift distortions are more difficult to correct); a more precise comparison depends on the assumption that the bias 
between galaxies and dark matter is scale-independent, and on other variables which are not of interest here, like h. 

The matter fluctuation variance in 8Mpc//i cells ag for a Q.m = 0.3 universe should be around unity to fit the cluster 
abundance (White, Efstathiou & Frenk 1993, Viana & Liddle 1996, Girardi et al. 1998). Wang & Steinhardt (1998) find, for a 
constant-w model, a general expression for as, that corresponds to as G (0.85, 1.25) at the 95% c.L, adopting our cosmological 
parameters. It is found that this is satisfied by < 0.1 , so that this can be taken as the upper limit on |/3| (see Fig. 9). As 
found analytically, the suppression of erg with respect to COBE-normalized standard CDM is due to the growth suppression 
in MDE. Another factor is that now the COBE normalization includes the integrated Sachs- Wolfe effect, no longer negligible. 
The rise in the CMB spectrum that we observed for small f3 induces via the COBE normalization a similar small rise in erg 
for the same values, as can be seen in Fig. 9. 

A fit to crg(/3) gives 

ag(/3) = ag(o)10(««" = -(«-^«^ (31) 

where crg(o) refers to uncoupled quintessence, and contains all the dependence on the other cosmological parameters, as well as 
on the exact COBE normalization scheme (I used here the Bunn & White (1997) normalization implemented in the original 
CMBFAST code). 



5 CONCLUSION 

Soon cosmology will benefit of high-precision data that will allow unprecedented accuracy in testing fundamental theories. 
Quintessence models add to the battery of cosmological parameters at least two entries, one describing the potential of the 
field and another its coupling to the rest of the world. So far, the coupling was arbitrarily put to zero, although we know no 
symmetry condition that implies so. In this paper we let the coupling be non zero, and investigated systematically all the 
possible trajectories in a flat space that lead to accelerated expansion at the present with = 0.7. 

The results for the homogeneous theory are applicable also to all Brans-Dicke models with a power law potential, since 
there is a direct correspondence between our constants f} and fi and the Brans-Dicke couphng constant ^ and the potential 
exponent n; from Eqs. (^, ^,|l^ we have in fact 

= ^{n-4)\. (32) 

For the fluctuations, the transformation that brings one from the fluctuation quantities in the Jordan frame ( the frame in 
which the field is coupled to gravity) to those in the Einstein frame (in which the field is coupled to the matter) is more 
complicated, and a complete treatment is still to be published. However, in the limit in which the fluctuations of the scalar 
field are small with respect to the fluctuations in the other components the fluctuation fields are conformally invariant; since 
the CQ field is almost homogeneous due to its light mass, this condition is verified for most of the universe history. As a 
consequence, it is conjectured that also the perturbative CQ results apply to Brans-Dicke models. The verification of this 
conjecture is left to future work. 

The main feature of the CQ model is the existence of a phase intermediate between the radiation era and the accelerated 
era, that we denoted tj) MDE. During this era the fluctuations grow less than in an uncoupled model. The ^MDE has three 
effects on the CMB: the spectrum at low multipoles is tilted, due to the ISW effect; the acoustic peaks are shifted to higher 
multipoles, due to the change in the sound horizon; and their amplitude is changed in a non-trivial way. On the power spectrum 
at the present, the main effect is the reduction of erg for large couplings and a very minor enhancement for small coupling. 

We found that the potential slope is not efficiently constrained by observations, essentially because the ^MDE is inde- 
pendent of /i. The coupling [3 is on the contrary constrained already by the present data, and is expected to be much more 
so in the near future, by at least an order of magnitude. From CMB and ag data we can derive the bound 

1/3] < 0.1 (33) 

which is stronger than the dark matter constraint of Damour et al. (1990) and not very far from the universal constraint of 
Wetterich (1995). 
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Figure 2. CQ phase space for values that Ue in the a region, fi = 0.1 and f3 as indicated. Tlic attractor a is the same as in the uncoupled 
model, but for /3 7^ there is a saddle for a non-zero value of the scalar field density. The phase spaces for the values of /3 investigated in 
this paper are qualitatively similar to the /3 = 0.5 case. The trajectory that falls almost vertically from top is similar to the background 
solution effectively employed. 
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Figure 3. Behavior of Qi/ (dotted line), Qft (dashed Hnc) and (thick Hnc) as a function of log(a) for ii = 0.1 and f3 as indicated. 
Notice that for CQ there is the transient regime 0MDE in which both the matter and the scalar field energy density are non-vanishing. 
Notice also that in this case, and for all values of /3 7^ 0, the matter-radiation equivalence occurs earlier than in the uncoupled model. 
For the small values of /3 used for the perturbation calculations, however, this is a small effect. 
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P=4.02, ^1=2.68 




Figure 4. Phase space of CQ for values that he in the bjvf region. There is a saddle cjjM at {x,y,z) ~ (0.12,0,0.95) that attracts the 
trajectory that falls from top, similar to the one adopted in the perturbation calculation. 
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Figure 5. Behavior of Qm (dotted line), Qr (daslicd line) and (thick line) as a function of log(a) for the same parameters as in 
Fig. 4 (label = 0.4) and for /3 = /i = 2.37 (label Wgfi = 0.5). Notice the transient epoch in which radiation and field share the total 
energy density (saddle cum)- 
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Figure 6. Evolution of the matter density contrast 5/a for three wavelengths, 27r/fe = A =1000 Mpc//i , 100 Mpc//i and 10 M.pc/h, for 
/* = 0.1 and various values of /3. Notice the growth suppression for /3 ^ 0. 
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Figure 8. Power spectrum A^(A;) = k^P{k)/{2Tv^) for the same solutions of type a of Fig. 7. The real galaxy data are from Peacock 
& Dodd (1994). 
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Figure 9. The Vciriance as versus p. The shaded area is the 95% c.l. region that matches the cluster abundance (Wang & Steinhardt 
1998). 
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